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OUTLINE OP RESEARCH ON OSCILLATING BOUNDARY LAYERS 


J . Coustelx 

National Office of Aerospace Study and Research 
(ONERA)i CERT, Toulouse, Prance 


1. Introduction 


Very diverse motives associated with the subjects of research /1^~ 

and application in fields, such as aerodynamics, hydrodynamics, bio- 
mechanics, meteorology and oceanography, have stimulated study of un- 
steady flows. 

We are particularly interested in the behavior of turbulent unsteady 
boundary layers in aeronautics, especially in helicopter and jet engine 
problems. The turbulence parameter is of fundamental interest in this 
limited field of unsteady aerodynamics. In fact, it appears that the 
analysis of unsteady flows can contribute interesting data for model- 
ing and calculation of turbulent flows. 

The most common unsteady flows are periodic and, quite naturally, 
they are associated with the basic case of oscillatory flow. We review 
recent work on this question, limited to that on the boundary layer. 

Other fundamental subjects, such as wakes, jets and flows in ducts, are 
being studied in several laboratories and also are very important, but 
they are set aside from this review. 

Except for the work ofKarlsson, it can be said that the development 
of calculation methods preceded that of experiments. This is partially 
explained by the necessity of having available means of acquisition of 
numerical data, while most of the means used so far have been analogic, 
this need being connected with the need of special statistical analysis 
of experimental data. Several experiments actually have been conducted, 
others are in progress, and they permit specification of the description 
of unsteady effects. They also permit a little broader discussion of 
the validity of available calculation methods. 


^Numbers in the margin indicate pagination in the foreign text. 
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Before describing these studies, we think it useful to recall some 
results on laminar boundary layers, which assist in understanding some 
aspects of the unsteadiness of the flow. 

2 . Laminar Plow 


Khile the equations of unsteady laminary boundary layers are known 
with great certainty, the understanding and analysis of seme questions 
have only recently been elucidated, in particular, by the use of asymp- 
totic expansion methods. 

The best known flows are oscillating flows, in which the average 
velocity can be zero or nonzero. Before describing the main results for 
these flows, it is advisable to stress some precautions to be taken, 
before using the classical boundary layer equations. 

2.1. Boundary Layer Equations (Incompressible) 

The Navier-Stokes equations quite clearly are the fundamental 
equations. In everything which follows, we are intest ed mainly in in- 
compressible flows, and this restriction implies conditions [Llghthlll, 
^3, 44] . In Isothermal flow, the two principal conditions are: 1. Mach 

number a little less than one; 2. the acoustical wavelength associated 
with characteristic frequency w should be large, with respect to char- 
acteristic dimension L of the system. With v as the characteristic 
velocity and c, the speed of sound, these conditions are described by 

v/c « ujL [c. A 

In this case, the fundamental equations are: 
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Discussion of the validity of the classical boundary layer equa- ^ 
tions, namely, 
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often brings numerous parameters into play in the unsteady case. V/e 
illustrate this with the simple case of a body of characteristic dimen- 
sion L, located in a fluid oscillating at frequency m. Uq is a typical 
magnitude of the amplitude of the velocity variation, and we shall see 
that the order of magnitude of the thickness of the viscous layer is 
(v/co)^^^* It uses a system of axes connected to the wall, and it in- 
troduces the folic wing dimensionless quantities: 

T-C06 Xiw/h 
^ z Y = t(coLVn)f/V^ 

Disregarding all effects of curvature, the two dimensional Navler- 
Stokes equations become 
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It then is clear that the solution of the problem brings into play 
two parameters: the Strouhal number oiL/U_ and the frequency parameter 

? Up 

coL /v . In particular, it is noted that, if mL /v is not large enough, 

p 

9 U term (Eq. 7) can be nonnegligible, and the classical boundary 
layer equations are no longer applicable. 


3 


In Gomo caGGG, only longitudinal length L 1 g not GUfficient to 
characterise the problem. For example, if the flow iG perturbed by 
a progresGive wave, the phaae velocity of which iG Q, it iG eeen that 
diGcuBGion of the validity of the boundary layer equationQ brings in the 
ratio of the wavelength Q/w to the thickness <S of the boundary layer. 

Therefore, by means of these examples, it is seen that the Reynolds 
number is not, as it often is the steady state case, the determining 
parameter, and that each problem requires detailed determination of the 
orders of magnitude of the various terms which appear in the equations. 

2.2. Oscillation with Zero Average Flow 

2.2.1. Stokes Solution 


The simplest case of an unsteady boundary layer is that of a 
parallel flow (v = 0) oscillating around an average zero. Besides, it 
concerns an exact solution of the Navier-Stokes equations, which is 
r>educed simply to 

H- -v> (9) 

It “ -bb " 

With the conditions at the limits, 

U - 0 y - n 

The solution found, as a function of reduced variable y=y(w/2v) 
is written 


with 


- cos Ulfc ODS Cu)t.^) 
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The roduoed amplitude Phase (j) profiles are presented in 

Pig. 1. An excess of amplitude in the boundary layer, with respect to 
the exter.^or, is noted. A large phase shift also is noted near the wall. 
At y » 0, the limiting velocity, i.e., also friction, is before it/h. 


u, » Au.cosut 

parallel flow 
61 6t 6y* 


T1 



U c AU COS(ul « if ) 


M « [(e-'sinTi)*-(1-e‘'cosii)’]''^ 
&u. 


lgs> a e-*slnn 
Ue*"cosn 



Pig. 1. Stokes problem. 


Practically, it can be /3 

considered that the thickness 
of the Stokes layer, is on the 
order of 5 ( 2v/o) ) • When the 

frequency increases, the unsteady 
effects due to viscosity are con- 
fined to a thinner and thinner 
layer. 

This problem can be consider- 
ed in a slightly different manner, 
by finding the perturbation pro- 
duced by an infinite oscillating 
plate in a flow at rest. Then, 
the conditions at the limit are 

V- 03.1 Uit cjs 0 

UA. O — 'M (SO 

and the solution is written 

— - c tos(ojt^u) 

0 


Thus, the oscillation of the plate creates a wave in the flow which 

1/2 

vanishes at infinity, the phase velocity of which is (2&)v.) 


2.2.2. Secondary Plow Induced by Oscillating Movement (Steady 
Streaming) 


When a body oscillates around a fixed position in a fluid at rest, 
it creates a secondary, stationary flow. 

Stuart has studied this problem in a reference system connected to 
the wall. The external velocity is in the form 
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The oheaa:» layer which develops at the wall Is analogous to a Stokes 
layer, and It can be studied by means of the following reduced variables 

U=U^V ccz'fL 

Is the current function) . 

2 

If (oL /V Is large, the solution Is calculated by means of the 
classical boundary layer equations. 


The dimensionless equations bring out the Strouhal number S=wL/v 
and. If this number Is very large, the solution Is found, by means of 
an expansion, as a function of the powers of Its Inverse a=S“ [43, 44] : 


( 11 ) 


In the 0 order, the solution Is classical. It Is that of Stokes: 


Ih, 


a 


(12) 


In order a, the solution Is more noteworthy, for stationary 
component Xg is not zero : 


(13) 


This secondary stationary movement Is Induced by nonlinear effects, 
due to the convection terms of the equation of motion. 


In order a, these terms bring out quadratic quantities, such as 

2 

sin wt, the time average of which is not zero and which necessarily gives 
rise to a stationary velocity. 
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It Is noted that the velocity induced at the boundary of the Stokes 
layer (^«>.<») Id nonzero. It is 


Ut s - J- Ug 

(iuJ "d'A 


(13') 


Thusj the problem is to determine how the average velocity becomes 

zero on the exterior of the Stokes layer. Two cases are distinguished, 

according to the order of magnitude of the Reynolds number R„»Un^/wv (this 

s u 

Reynolds number is made up of the order of magnitude of given by 
formula 13', U^'X/Uq^AjL and of the characteristic length L) . 

If Rg>>l (Pig. 2), a boundary layer, in the very classical sense 
of the term, develops below the Stokes layer, and its thickness fig is 
on the order of L / . The ratio of 6^ to thickness 

6g'\y(v/w) ^ of the Stokes layer is on the order of LUnOvg/Cv/w) 

= ojL/U «S. It is very much greater than 1. If R„ ^^“is on the order 
of 1 or less than 1, the external layer is of a different nature. It 
should be studied, by using the linear Stokes equations of slow move- 
ment . 





Pig. 2. Oscillating circular 
cylinder. 

much less, it preserves its individual 
of these two layers is 


2.3. Oscillation with ?Jonzero 
Average Plow. 

When the boundary layer 
develops in the presence of an 
external flow oscillating around 
an average nonzero Uq, two solu- 
tions play a fundamental part: 
the Blaslus solution and the Stokes 
solution (Pig. 3)* If the .thick- 
nesses of these two layers are of 
the same order of magnitude, the 
Stokes layer is fused into the 
stationary layer, but, if the 
thickness of the Stokes layer is 
nature. Now, the thickness ratio 
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-V. . S-"* 

Tho order of magnitude of the Otrouhal number, therefore, io a funda- 
mental parameter in definition of tho otructure of the boundary layer 
in ouoh a case. 


y U» j» Uc. ‘Ut, 5« ut 


Blastus S , JLi 


K_ * ptokoo layW 

it./iilV" ^ 

6 iu,/ 

Fig, 3 . Oscillating flow on flat 
plate (laminar). 


The conditions at the limits are: 


The analytical oolutlono of ^ 
this problem generally have been 
found, on the assumption that tho 
external velocity fluctuation is 
small 


^•s b. (/I + i, with5 «'^ 


(14) 


S'VmO U=0 •5C'>Q 


•4X, U- 


= Oe 


nf, a O > 


The solution was found [Llghthlll, 373 , by expanding the velocity 
by powers of e . By introducing the reduced variables 


A/ 


and the reduced current function 


(uj/ vtp* ) ijj 


the first order solution has the form 




(15) 


If the Reynolds number Uq x/v is large enough, solution P is the 
Blasius function P., a function^ of n=U„ alone. 


In order e, solution f can be determined numerically, whatever the 
reduced frequency 5b [^5-48] . However, analytical solutions are available 




'V 

f©r extromo valuoo of x. 


V/hon X io small, f io expanded by bho powers of x. The colution 
hao the form 

H MP '1- <f (10) 

whore fj.^ lo the solution oi^ a differential equation whleh uses 


For n“0, a quaoistatlonary boundary solution is obtained (x‘>-e) . 
fQ is obtained by writing that Blasiuc function F^ lo the solution at 
each Instant, The utilisation of expression (1^0 for U_ and expansion 

w 

by a small perturbation results In 


therefore, 



(17) 


When the Strouhal number is large, the boundary layer has a two 
layer structure. Near the wall, a shear layer of the Stokes layer type 
develops. Its thickness Is on the order of (v/w)^'^^. This Internal 
solution Is a function of y»y(co/v)^'^^, and It Is written 


The principal term is the Stokes solution. 

The external layer has a thickness on the order of x/(U x/v)^*^^, 

v» 

and Its solution Is found, as a function of 2- "k" ' 

It is written 

|s ^ .(4.L)/(5^)''' .i(F .^1 F')//i * ) 

Pigs. 4 and 5 give some summary results. Evolution of the depth 
of displacement and wall friction and solutions for Intermediate values 
of the Strouhal number have been obtained numerically [45-^73 . In 
particular. It is noted that the phase of the coefficient of friction 
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tondo tev/ard while x/lJ„ 
fcendo fcov/ardo Infinity, ao 
in fch*^ 3t©Heo oolution. It alee 
Id noted that, when the Gtrouhal 
iD largo, tho extornal region of 
bho boundary layor io cubjeobod 
to puloation ao a unit, v;hon 
tho StokOD layer io very thin. 
Solution (19), in fact, ohowD 
that, in thio area, there io 

t girtust 

Consequently, product 

OO 

remains conotant 
over time, for the very small 
region in tho vicinity of the 
wall contributes little to 
This explains the behavior of 
tho amplitude and phase of 6^ for large values of wx/U^ . 

G 

2 

In order e , a stationary term Un appears in the solution. If /5 

, 2 k 1/2 2 2 " 

numbers R„=(vw/Uq ) and e U^L/v are close to one, the boundary 

layer equations are applicable to calculation of C^5-^73 . With 
available experimental and numerical results taken into account, it 
appears that these secondary velocities are small, even for large ampli- 
tude external volocltleo. If this so, an important simplification re- 
sults, since the average profile is the classical solution of the sta- 
tionary equations. Sometimes, with an average pressure gradient, espe- 
cially in the vicinity of the separation, the secondary velocities can 
become Important. 

2.H. Eff e ct of Phase Velocity 

; The response of the boundary layer to a perturbation in the form 

slo l/J ( i m ( 20 ) 

^ 10 


UesUo(Ucsinwt) 



Pig, tu Oscillating flow on flat 
plate (laminar) j amplitude and phase 
of <S^ (from [293 ). 



Fig. P* Ooeillatins flow on flat 
plate (laminar) j amplitude and 
phaoG of wall friction (from [29, 
45-'l73). 


deptndc greatly on the cpeecl of 
eonveotien &f the wave Q. M.II, 

Patel har deno an oxporimonfcal 
ctudy, by producing a oyotem of 
odddoD ooeaplne fcon a owinginB 
ohutter (Fig. d). Loeatod in the 
flow at voloeity U^, the oadiec 
move by eonvootion at velooity Q» 

0.7? Ug, 

In the first order, the rosponoe 
of the boundary layer to the por- 
turbatiOii is in the form 

U 5 bo Sir^f (21) 

The experimental rooults indi- 
cate that the profile of continuous 
component is not affected 

by the oscillation. It remains a 
Blaoiuo profile (U„ is Indepondont 
of x) . Actually, oven if there were 
nonlinear effects, they would be 
very weak, for the reduced amplitude 
is 10/^ at most. 









.cd: 

'-ai^ 

1 


miJf* 


The Ut/Ut amplitude and phase 
(j) profiles have some peculiarities, 
especially around Strouhal numbers 
wx/Uq in the vicinity of 1 . The 
amplitude has values very clearly 
higher than 1, on the order of 1.6 
(Pig. 7), while, in the case de- 
scribed in the preceding section 
(which corresponds to Q»-«>), the 
amplitude becomes not so much more than 1, around 1.1. For some 
vaD”ss of o3x/U^ , the amplitude profile, moreover, has values clearly 


rx--_':irr 


Pig. 6. Experimental assembly 
used by M.H. PatcTi the laminar 
boundary layer was studied on 
auxiliary plate a. 
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less than 1 toward the outside of the boundary layer. The theoretical 
elements developed by Patel, which are an extension of the methods of 
Lighthill, are valid for extreme values of the Strouhal number (wx/Uo ->-0 
and cox/Uo^-><» ) , do not allow explanation of these peculiarities, 'ccord- 
ing to Patel, they must be associated with the fact that, when wx/Uq^= 
0(1), the inertial and frequency terms are of the same order of magnitude. 


U,a SinO(t-Ji) 

Q' 

U s Uj ♦U, sin[w(t-^)*4>] 


Q = 0.77Uo 



1. lov: frequency sol. 

2, frequency sol. 


Pig. 7- Laminar boundary layer 
of flat plate perturbed by a pro- 
gressive wave (from [33-36]). 


The phase profiles Indicate 
that the velocity in the boundary 
layer always is slower than the 
external velocity. 


Two parameters compete to 
determine the phase distribution: 
velocity Ug and pressure gradient 
-1 3P. Now, the expression for 

p 

this gradient is 


if. - 

f &0. ■ Ifc 


IS tWt) 

_.y SirU^to/^. 
a Q ' Q/ 


( 22 ) 


At low Strouhal numbers, the ^ 

velocity parameter prevails. The behavior of the boundary layer tends 

to become quasistationary, and the phase shifts are small. When 

0 ) Increases, the pressure gradient parameter becomes large (Eq. 22). 

If Ut /Un is small enough, the first term on the right side of (22) is 
■^e e 

dominant and, when convection velocity Q is smaller than U^, the pressure 
gradient term lags 90° behind . Besides, the interior of the boundary 
layer is more sensitive to the pressure gradient than the exterior, 
since the inertia is less. Therefore, this explains the negative phase 
shift values. This same reasoning explains why It Is not paradoxical 
that the phase shift Is positive when Q-x^ (see Section 3.3). 


A singular situation results when Q=Uq, for the pressure gradient 

is zero in the first order. The solution proposed by Patel, for large 

values of oix/Uq , effectively Indicates a change In behavior of the 
•8 

amplitude and phase profiles around Q=U^ (Pig. 8) . 
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3. Turbulent Flow 


U,s Us+U,, slnw(t-X) 


3.1. Experimental Reoult:s 



In turbulent flow, under- 
standing of unsteady effects is 
much less advanced than in laminar 
flow. The experimental efforts have 
concerned some particular points. 
Several studies have been devoted to 
oscillating flow around an average 
zero, with the effect of wall rough- 
ness. These studies often have been 
motivated by shear wave problems in 
the depths of the sea [21, 22] . The 
effect of pulsation on a boundary 
layer of a flat plate also has been 
studied [9, 23^ 20, 11], but system- 
atic analysis of the various paramett 
which can affect the structure of the 
turbulence is still far from being 
achieved. In continuing study of 
laminar flow, M.H. Patel also has 
carried out an experiment in turbulent flow, on the effect of the phase 
velocity of the perturbation wave but, there also, the field of study 
remains largely open. Finally, very recently, several studies have been 
directed towards the effect of an average pressure gradient [3, 12-15, 
38, ^1]. 


Fig. 8. Effect of wave convec- 
tion velocity (from [33-36] . High 
frequency theory wx/U =5). 


3.1.1. Oscillation Around an Average Zero 

In Europe, studies of a turbulent boundary layer induced by oscil- 
lation imposed on the flow have evolved mainly in Denmark since the 
1960s. Recently, Jonsson published an article which focuses on the 
group of studies carried out at the Technical University of Denmark. 
There, he also presents some outside results. The aim of the analysis 
is to demonstrate the existence of universal relations for the velocity 
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profiles, and to establish practical rules, especially for the coef 
ficient of friction. The results are mainly for a rough wall. 


Accordingly, Johsson was interested in flows, for which the exterior 
velocity is in the form 

Ml = wb 

In the boundary layer, at the first harmonic, it is written 

U s. cos(ojb + 'f) 

‘ ( 24 ) 

We note that, since the flow is turbulent, it is advisable to 
define the average u correctl;^ . This velocity is the average of a 
set and, when the flow is periodic, it also can be the phase average. 
Therefore, u is an average of the instantaneous values recorded in 
different cycles, for the same phase angle wt within 2 tt . 

Jonsson studied a deficit velocity, which he defined in complex 

form: 



In a stationary turbulent boundary layer, conclusive progress has 
been made, due to the existence of a recovery region, where the velocity 
changes logarithmically. Under rough conditions, the rule of the wall 
is found by using the following reduced variables 

0 / and K 

where U is the friction velocity, U =U (Cf/2)^'^^, and k is the height 

IT* J/ s 

of the roughness. 

Jonsson attempted to verify such a logarithmic rule, by consider- ^ 
Ing velocity U and friction velocity at the same Instant. 

l4 


;s»rss“- 

The friction velocity was determined from measurement of the 
velocity profile and by means of the summary equation of the mo- 


mentum. This equation Is quite simple. It results from the local 
equation which. In the case of parallel flow, is written 

“ bt dy p C 26 ) 

and it takes the form Iz - u") 

e ■* 0- (27) 



I^lg- 9* Turbulent boundary layer oscillating 
around an average zero (from [21]). 
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Pig. 9 shows that, during a large part of the period, the measure- 
ments are grouped around the same curve as in stationary rough flow, of 
the equation U/Uj^=5.75 log(30y/k). Large deviations develop in the 
neighborhood of times when Ue is nullified. 

Determination of the eddy viscosity shows that, In these same 

periods, has negative values, because of the inertia of turbulent 

friction -p<u'v’> , with respect to the constraint Meanwhile, these 

Dy 

periods remain quite limited and, for practical calculations, the utili- 
zation of a hypothesis of eddy viscosity should not be rejected a 
priori. In any case, it must be noted that the frequency imposed on 
the flow remains low with respect to the estimate which can be made of 
a characteristic frequency of the turbulence. This can explain why some 
classical results of stationary flow remain valid. 

In a less convincing way, Jonsson also establishes the existence of 
a universal rule of deficit velocity, in the form 

(bis, by convention, the closest point of the wall, where the maximum 
velocity equals the maximum external velocity) . 

The existence of such a rule and of the rule of the wall, in the /_^ 

form U/u^=f^(y/k) , implies that of a rule of logarithmic recovery. 

Jonsson proposes 

?.?5 

The extent of the recovery zone diminishes, when the ratio of the 

amplitude of the movement to the outside a^ to the height of the 

roughness k decreases. This zone disappears, for values of a-i /k on the 

e 

order of 30 . 

It also must be noted that the existence of a rule of the wall 

U/Uj,= f^(y/k) implies that the phase shift of the velocity in this region 

remains constant. In fact, it is enough to allow that is in the form 

U =U cos(wt+bn) , and the following is obtained 
r rm ^0 ’ “ 

Cos(ojl;^^o) 
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Experiment showe that, in the recovery region, the phase shift 
remains quite constant, but it appears to evolve in the viscous region 
under the layer. This then denies the existence of a rigorously uni- 
versal rule in this underlying layer. 


In order to permit complete construction of the velocity amplitude 
and ph^ase profile, Jonsson found a representation of the phase profile 
of the deficit velocity He proposes 

Ti - 1 

Practical rules likewise have been proposed for the parietal coef- 
ficient of friction and the boundary layer thickness b. For thickness b, 
he proposes, under rough conditions. 

Two conditions are distinguished for the coefficient of friction 



< 4,6 


o.i'5 

«C! 0,\ 

A 






> -^.6 


0, 0 <305 

A 


h‘ 



Jonsson also has proposed formulas for the smooth condition 

: 0,03 Hc;'is£± 

fu‘ ' ■>< 

0 , - U 4 P 

■* 

3.1.2. Oscillation Around an Average Nonzero Constant 


In this section, we are Interested mainly in the case of a boundary 
layer of a flat plate, perturbed by a sinusoidal oscillation of the 
external flow 


u, - cot 


( 28 ) 
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At the first harmonic, the velocity in the boundary layer is In the 


form 

u s .t cos(ojhf) +... 

(29) 

As in Section 3.1.1. , U is the average of the set of instantaneous 
velocities . The difference between the instantaneous velocity and the 
average of the set is the turbulent fluctuation u' . The moments, also 
defined by the averages of the set, are designated <u' >, <u'“>, and 
the turbulent pressure is -p<u'v'>. 


This problem is the basic circumstance of numerous practical situa- 
tions, and many numerical solutions have been proposed over several 
years. For a long time, the only supporting experiments were those of 
Xarlsson. Recently, others have published [9^ 20, 11, 12] . The table 
below gives the principal characteristics of the various cases studied. 
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In this table, X designates the distance of the plate on which the 
boundary layer is studied to the leading edge or. In the experiments of 
Houdeville et al and Cousteix et al, the distance from the fictitious 
origin of the boundary layer, calculated from the experimental data, on 
the assumption that the average characteristics (in time) of the boundary 
layer obey the stationary laws of the flat plate. A characteristic 
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frequency of the turbulence has been defined by f^=u'/6 , where u* 1e 
a characteriDtic velocity of the turbulence, for example, /crr2 . 


In the experiments of Karlsson, the moasuroments provide average 

velocity profiles and profiles of in phase components U^^cos^ and out 

of phase components U^^sincj) . The experiments of Charnay and Melinand 

give profiles of the average velocity, the average of the set U at 

different instance of time in the period and some longitudinal and 

transverse turbulence intensity profiles. In the experiments we have 

performed, we measured the profiles of the average velocity U, which 

are analyzed by means of harmonic decompocition, and we reported the U^, 

Un and ^ profiles. Profiles of the longitudinal intensity of the 
2 1/2 

turbulence <u'^> ' and turbulent pressure -p<u'v’> also were reported. 
The density of probabilities of u' also was determined, as well as the 

h 0 0 

cblateness factor P=<u' >/<u > . 


Average Velocity Profiles 

The most evident conclusion of the experiments of Karlsson and 
our experiments is that the average ’velocity profile is practically 
unaffected by the unsteady nature of the flow. However, it seems that 
the experiments of Charnay and Melinand contradict these observations. 

They particularly indicate an amplitude effect. Actually, for a 

of 0.2 (f=l8.5 Hz), the average velocity profiles clearly are 
different ("squarer") than those obtained for Un /U^ <0.11. For this 

‘ J.Q Oq 

same value of the reduced amplitude, an 80^ increase over the station- 
ary case of the average thickness of the boundary layer has been measured, 
as well as a 13 % increase in the coefficient of parietal friction 
(X=0.9 m, f=l8.5 Hz). The observed effects on the average properties 
of the boundary layer seem all the more surprising, that the experimental 
conditions are approximately in the range of those studied by Karlsson. 

In any event, if the effect of a turbulent perturbation of the 
external flow [2, 8] and that of a harmonic perturbation are compared, 
as Charnay and Melinand note, the effectiveness of the latter on the 
average characteristics is very weak. The fundamental differences be- 
tween the two types of perturbation must be emphasized. In one case, the 


energy is distributed in a speotrum similar to that of turbulence, and 
it is a band in the other. We note that, in some cases studied by 
Karlooon, the frequency of the perturbation between the range of char- 
acterJotie frequencies of the turbulence and the effect on the average 
profiles is not so large. The tridimensional and unidirectional 
natures of these two types of perturbation also must be compared. 


Perturbation and Phase Ampl itu de 

Qualitatively, the amplitude profiles we obtained and those of 
Karlsson have the same trend (Pigs. 10, 11). In particular, they 
exceed the outer value in the external regdon of the boundary layer. 
These excesses can be on the order of 15 to 20^. V/hen the Strouhal 
number increases, a layer similar to the Stokes laminar layer develops, 
of less and less thickness in the vicinity of the wall, and in which the 
amplitude variations are confined. 
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Pig. 10. Turbulent boundary layer in oscillating flow. 




Fig. 11. Turbulent boundary layer 3n 
oBoilXatJng flow*. /U^ «3^*^?5 -wx/Uq= 

3.92-R ^=0.88-10^ (from [23]). 

©X 

The phase <}> profiles are not given explicitly in the article of 
Karlsson, but they can be determined from the in phase component 
U^cosd and the out of phase component U^sincj) profiles. With the scale 
on which these data are reproduced taken into account, some inaccuracy 
results from these analyses. Meanwhile, it appears that the results we 
obtained completely ccrfirm these interpretations. Generally, it is 
observed that the phase is positive in a zone near the wall. As in 
laminar flow, this can be explained by the fact that the pressure gradient 
leads the external velocity by 90°, and that the less inert boundary 
layer in the vicinity of the wall responds more rapidly than the external 
flow. Toward the outside, significant negative values have been measured 
for some Strouhal numbers: 2 . 95 , 3.92, 5.89 [23] j 4.6 [11, 12]. For the 
other Stx’ouhal numbers studied, the recorded values of c|) toward the out- 
side are generally within the experimental error. Similar observations 
have been made in laminar flow (measurements of Hill and Stenning, cal- 
culations of Parn, Arpaci and Clark) . As a function of the Strouhal /lO 
number, the lag of the phase maximum in the boundary layer is at a 
maximum for values of wx/Uog, on the order of 2. 5-3. 5. 

It also is noted that variations of c(> are established in a zone, 
which becomes thinner when the Strouhal number Increases. At large 
Strouhal numbers, the outside of the boundai'y layer pui.ses as a unit. 

The phase shift is very small, and the amplitude is very close to the 
external value . 
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i^napo Ehaft Relation Near the VI oil 


Pinall,y, it Ic noted that the phocc near the vmll han a mQjdmut« 
and dccroaooD at lower values of j;. It is noted that the existoneo of 
a logai’lthmie recovery resJon ImplJes that the phase shift must he 
oonstant and equal to it£ value at tie wall [^113 . It also can be shown 
heru, by ascumlng that thoro f.s a. universal rule near tho wall, in the 
form 


) 


(30) 


For this, we use the following complex not at lone 

. * 

y To 4 T. 

, Cwf 


<JL a W* ^ -Ua a 


vjith L, . Un . Y « IJ 

- - eg- Ig- To' 


o 


actual and and complex ► By assuming that 
1 compared b>': , y v.nd P , the rule ■u’' “ 


li, , Y't ar.U U. arc small .... . ... ...>. 

4.0, ' i. .( ^0 

f (y ) can be developed, and the following are obtained 




'0 ^ ^ y dy h I 

‘Ut,. rj(i‘ V I V 


(31) 

(32) 


Formula (32) indicates that the phase t|) of u is dndcpei^aei.t of 
y. It is such that 

y / !£i *. iO«.\ 


^ 3f = 


^ ^ r, d.J 


1/2 

Therefore, it is equal to the phase of YU„=(r^/p) and therefore, to 

© P 

that of the wall friction e , 

P 

This conclusion is compatible vcith experiment, otuly in the region 
where the phase passes its maximum, and it can be thought that the 
interpretation of the variation of (j) very (.'lose to the wall requires 
sharper analysis of the viscci.s underlying layer. 
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Mearwhilt , it er.i. bo noted that the maxirr.um of ^ qqqixvg In a hof.ton 
'./here the* logarithir.lc rule Ic ectabliDhod in ct at ionary f»©w fexeept for 
oxti'cnie very high v&Iugg of the Strci'hal mur.ber L*t'’d?ed by Kerlcoon, 

S»77»07). V/e grant that the nittlonai'y rule romaino valid at 
oaeh Inetant 

U"*" - i'll 

4 0 f..ov 

On th.« sUj *t hand, DimiHarity eolutirinr; hevi beet extended to the 
ur steady ease, asBumlng a dv.flcH vcv'.eoity rule cf the form (tj^-U)/Uj,® 

P' (y/ 6 (x,t) ) , and by aeoepting a quasictatlonary mixing length cchcme 
Dil. It then is shc\in that friction can be calculated at each instant, 
by the relationship 

with 
and 


(S s iili H 2 it- 


(3^0 


A small perturbatiot hypethesis is then made, end it is written, 
in complex rotation 

a a Oo + 

Ofc = tec, » Ue 

r Z t * I eM 

Su s S,„ + 

H s Ho + K, 

0 - 6. + 


The expansion of relationships (33) and (34) then leads to 




+ (e-»-''J>)A;;^05'f^^-3)/'0Cos'f£> 


with 


>- 

It 

(2 . w \ 
1 c'' / 

‘ M, 

C: 

V. 

U 

5 = 

(•* - 

|.± 

3: 

B. 


' ACo'«y 

' Ho 

J 


(35) 


23 


c.i* 


(p 



gj ami ip^ avc plKijtr cl ohU Q, tml Iq 
ara 0, recJuaccl by tin* an-pJ*Ut4t' of 0^,"' 



I Sjil /Sf0 

^•>c f 



!ul /Oo 

»»awKU» iwii.jmjii 

fJH 


Lr: Cm amplifcuSoo 


Ti'fC right hanj term ©f Eq, (35) aari ho detoi^mlnctl from oxporJmcr tal 
data, fho value whleh. the* veleeitSr phrce rheuld hovi. 3rt the logarlthmie 
rogieK then ic actormSrecI. Thio phaee cheulh correorxnd te the maxlmurc 
value detc.rii’ined experimentally, v/hioh alcr.e io ceropatiblo with 
the theorotleril reculfc of $*»ooto. 


Cemparleeri of the theoretical v&luee determined v;ith the 

experimental, values is pi‘er.ent<-'d in rig. IT. The rooults cf Karif.scn, 
ur.fortunoiely . oould not be used, teeause the Inform&tion ic lnouffl~ 
dent to calculate the right hand term of (25). Our results have been 
plotted in Pig. 12 [20, 11, 12] , Likewiso, those obtained v?ith an 
average} pressure gradient [12-153 and those of M.H. I’atel, who has 
studied the effect of th^.- convection velocity of the perturbation wave, 
also have been plotted. Thu resulting agreement is altogether good, ever, 
for relatively high values of the Btrcuhal number, on the order of 6 to 7. 



Effect on Turbulence Charaettr- 

1st iCG 

In the exporir.erts of 
Karisson, only the average time 
value of the Intel si ty cf the 
turbulence was measured. It 
is impossible to detect a pos- 
sible effect of frequency on 
these profiles, espeolaliy with 
the experimental error and the 
fact that higher order har- 
monics then i .are Included in 
the measurement taken into 
account. 


Fig. i2. Ve.i'l fief t-on of propesgd, re- 
lationship for phase shift in the 

semiJ. ogarithniic region. e 

2H 


'ihv vjyik of Chawicy &r>cl P,c3-»r*mia i«velvcd rtwly t** ♦ft. iTt:f Lo»;na*" 
.'•sy f»nC the- voloeftico en the inaiile tiM cmtaicic of the turbulont burotc. 
It le ev’idoht that, foi* the ctst.torifiry t*r-at , the Indentationn ©f the 
boundary are deeper end that, ptebr.b3j, the ortgln of the buret© ic deep- 
er ineide the boiimlery iayi. r. Or-fortunGtoly , thic etudy Id of toe 
eutnmary a nature, for ana^yclD ao a funotien ©f the phaoo hao net been 
oarried out end, w dor theoo eoncUticnc, it In dH’fieult to make a 
dlctJnct-ion botwoen the foreod variation© imporod on the free In-undary 
by the pulcation of the flov; and its random vaidotion, v?Pieh probably 
Id a function of the phaoe. 


.lore aotriSlid mcar (,* eimitD of the buhavlor of certain turbulent 
quantltiec have been reported by Gouotoix ot al, mainly in the cecond 
aoriflruratien, which correopondo to the hlghoet Strouhal numbfuv, 
Moauurementa of thQ oorrelatlon aocfflciont xv *^>) , 

which have an important role In modeling turbulence, have boon carried 
out. They ohow that thle coefficient ev elves in the came way as in a 
stationary bounda.ry layer vPig. 13). Hear the wall, within the moasure- 
mont cin'or, it can be considered that it is oonotont, with a value on 
the order cf 0.45-C.f. Turbulent pressure measuromente also have per- 
mitted detormdnaticn of the; experimental development of the mixing 
length il==( -<u' V>) This operratlcr is c tit critical, for it 

requires calculation oPthc; derivative 3U/3y. however, it shows that, 
in the entire period, the mixing S.ength evolves in the boundary layer 
in a very classical manner, except In a sene located In the vicir.itj c..t 
the external velocity maximum, where large variations are noted (Pig. 14) 


It also is of interest to analyze the evolution of the oblateness 
factor P«<u' ^>/<u' especially if it is Interpreted as discontinuity 
factor Yj t>y using the relationship y^’S/P- b’e recall that y, in unsteady 
flow, is a functicn of uhe distribution cf the free boundary of the 
boundary layer. Tho yiy) curve (Pig. 13) depends on phase, and this 
dependence esacntlaliy expregses the forced variation cf the boundary 
layer thiePness. Actually, if y is plotted at yy-Q^^ (the value of y 
for v/hich y“0*^’)j the points are grouped, and they form a curve sim- 
ilar to that given by Klebanoff, in the c&Be of the stationary boundary 
layer of a flat plate, This tends to shew that the purely random 
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Fig. 33 . Puloed turbulent boundary 
layer, correlation coefficient. 




Pig. l4 . Pulsed turbulent boundary 
layer, mixing length distribution. 


variation of the free boundary 
in not affected by the unsteady 
nature of the flow. 

Other measurements reported 
in detail by others [11, 12‘J 
tend to indicate a quasi- 
stationary behavior of the tur- 
bulence. However, it must be 
noted that the measurements were 
made under conditions, such that 
the frequency imposed is quite 
clearly less than the character- 
istic frequency of the turbulence 


An estimate can be made of 
the conditions under which these 
two frequencies are of the same 
order of m.fgn-' t ade . A charac- 
teristic frequency of the tur- 
bulence is f;p=u*/(Si where u' is 
a characteristic velocity of 
the turbulence, for example, 

(u ' ^) It can. be estimated 

tha t the order of magnitude of 
u'2 is given by (rp/p)/0.3 
(it is assumed that, near the 
wall, -u'v'-O.Su'^) . By using 
the flat plate formulas 0^-,= 
0.0368/R^/-/^ and 6/,x=^C 
it is found that f,p/fc;_ 2 'n'R^l/ 2 ^ g , 
f,p and f are of the same order 
of magnitude, if S=2TrR.^^'^^^ . 

For Reynolds numbers of 10° to 
7 

10 , this condition is satlsfjeo 
for very high Strouhal numbers, 
on the order cf 20. Under such 




Fig. 15' Pulsed turbulent bound- 
arj' l;;yer, ob lateness factor. 


conditiono, unsteady effects develop 
mainly in the underlying layer, and 
the bursting phenoiticnen can be affected. 
For higher frequcnciec., It can be 
asked whether a frequency does not 
exist, beyond vdiich the unsteady ef- 
fects on turbulence become very v;eak 
since, then, the unsteady layer is 
limited to an almost purely laminar 
zone near the wall. 

3.1.3* E ff ect of Wave Convection 

Velocity 

M.H. Patel C3^-h-36] has developed, 
in turbulent flow, a study similar to 
that described for laminar flow in 
Section 2.4. The external velocity 
is in the form 

with q _ o,?7Uo^ 


The general conditions of the experiment are summarized in the 
table below. 



It is noted that the imposed frequency is quite clearly different 
from the characteristic frequency of the turbulence. 


The observed effects on the behavior of the average velocity, 
amplitude and phase profiles are qualitatively the same as in laminar 
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flow. The average velocity profile is unchanged from the stationary 
case. The amplitude exceeds its exterior value (Fig. 17), and the re- 
duced amplitude maximum /U^ increases with the Strouhal number 
(Pig. 16) (it is 1.77 for ®S=6?65) . 
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Boundary layer of flat layer flow (Pig. 16). 


It is noted that the 
Strouhal numbers studied are 
larger than those in the study 
conducted of laminar flow, and 
that the singular results ob- 
tained around S=1 were not ob- 
served in turbulent flow. 

As in laminar flowy the 
phase shift in the boundary layer 
always is negative (Pig. 17). 

The absolute values of the phases 
sometimes are smaller in turbulent 

As was discussed 


perturbed to progressive wave (from section 2.1t, two parameters 

come into play for the phase 

shift behavior: velocity and pressure gradient. Given that the tur- 

bulent boundary layer is less sensitive to the pressure gradient than 
a laminated boundary layer, it is logical that the phase shifts are 
smaller than in a turbulent boundary layer. 


Finally, we note that, in a large part of the boundary layer, <j> 

remains nearly constant (Pig. 16) . This result is completely com- 

"t" + 

patlble with the existence of a universal rule u =f(y ) (see Section 
3.1.2). We recall that the values of these phases were compared to the 
theorelrical values, calculated by Eq. (35) (Pig. 12). 


Likewise, M.H. Patel measured the longitudinal Intensity of the 
average turbulence over time. No effect of the pulsation was recorded. /13 
Within experimental error, the profiles are Identical to those of 
stationary flow. 
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Pig. 17. Turbulent boundary layer of 
flat plate perturbed by progressive 
wave (fromM.H. Patel). 
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3.1.^. Effect of Average 
Unfavorable pressure grad- 
ient 


Several experiments now 
are under way or planned, to 
study the effect of pulsations 
on a boundary layer exposed to 
an average pressure gradient 
[ 38 , 41 , 22 , 3 , 12-15] . The 
principal characteristics of 
these experiments are summarized 
in the table below 
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In the above table, X«0 corresponds either to the leading edge of 
the plate on which the boundary layer was studied [2^, ^1], or to the 
fictitious origin of the turbulent boundary layer, determined from the 
experimental results at the first measurement station [38, 12-15], or 
estimated [3] . 

It Is noted that the use of a hydraulic cell [3] Is a good means of 
easily obtaining high Strouhal numbers and pulsation frequencies, which 
are In the range of the characteristic frequencies of the turbulence. 

In the experiments reported by Kenlson, no unsteady effect was 
noted on the evolution of the average characteristics of the boundary 
layer over time: velocity profiles i coefficient of friction^ density of 

momentumj shape parameter. Intensity of turbulence even In vicinity of 
point C^=0 (Pig. 18) . Upstream from C^=0, a zone develops, where the 
velocity periodically Is negative. This zone remains thin, and the 
rough thickening of the boundary layer Is observed practically In the 
same zone as In stationary movement. 

The amplitude and phase profiles also were measured by Kenlson. 

The phase profiles Indicate tendencies associated with phase shifting of 
the pressure gradient ^ with respect to the external velocity. The 
phase shifts In the boundary layer are more pronounced than in the case 
UQ=cste. This can be explained by the fact that. In the vicinity of the /l4 
breakaway, the velocity near the wall Is much lower, and the inertia, 
therefore. Is much weaker. 

The experiments of Schachenmann and Rockwell In a conical diffuser 
with a 6° total aperture Indicate that the average characteristics of 
the boundary layer are not affected by variation of the pulse frequency. 

The average pressure recovery coefficient Is not affected either. Pre- 
liminary tests conducted by Binder et al with much higher Strouhal 
numbers lead to identical conclusions, as long as the boundary layer is 
not separated (2a<7°)* Por larger apertures, very noticeable Improve- 
ments in yield are obtained by Increasing the frequency. It must be 
noted that. In these tests, the external velocity frequency and ampli- 
tude varied at the same time, and It Is impossible to separate the 
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effeota due to frequency and 
those due to amplitude, Very 
strong relative amplitudes 
were obtained up to 1.24. In 
addition, this problem brings 
Into play coupling between the 
external flow and the boundary 
layer . 

The tests reported by 
Simpson were carried out at low 
Strouhal numbers and at low 
frequencies, compared with the 
characteristic frequencies of 
the turbulence. The first re- 
sults published In this study 
mainly concern the average set 
velocity, amplitude and phase 
profiles. At the station pre- 
sented, it was found, in par- 
ticular that the average set 
velocity profiles all have a 
semilogarithmlc region, in 
which the phase Is practically constant. The great thickness of the 
boundary layer has permitted precise determination of the phase varia- 
tions in the underlying layer. These soundings confirm the very high 
sensitivity to unsteady effects, and they show that phase <j) varies very 
rapidly . 

The first results we obtained on the effect of unsteadiness on a 
boundary layer subject to a pressure gradient were P^r’esented in detail 
during this colloquium. In particular, data obtained by means of a 
laser anemometer in a zone with back flow are described there. As this 
zone develops j the unsteady effects become larger. For example, during 
one period, the shape factor at the last measurement station varies 
between 2 and 5- When the shape factor is raised, very noticeable modi- 
fications of oblateness factor' P are noted. In particular, P Increases 
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perturbed by progressive wave (from 
Kenison) . 
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in the vicinity of the wall . Meanwhile, it wao verified that the 
average set velocity profiles at each instant of the period can be rep- 
resented by profiles, drawn from a family established in the stationary 
state, in the outer region of the boundary layer (Pig. 19). This was 
observed, as long the Instantaneous profile did not have a negative 
velocity (H<2-2.3)i up to Strouhal numbers on the order of 6 to 7. 

Likewise, the existence of a semllogarithmlc region could be observed 
at each Instant during the period, the extent of which decreased, as in 
stationary flow, when the shape factor increased. It practically dis- 
appeared when H was on the order of 2 to 2. 3* In this region, the 
velocity phase angle in the boundary layer, with respect to the external 
velocity, has an extreme, the value of which has been compared to the 
theoretical value given by Eq. (35) (Pig. 12). Good agreement is obtained, 
up to the time when the average velocity of the set becomes negative in 
a small portion of the period. Here again, rapid variations of the 
phase angle were recorded near the wall. 



Pig. 19. Pulse turbulent boundary 
layer with average pressure gradient ^ 
comparison of Instantaneous velocity 
profiles to theoretical family. 


3.2. Calculation of Boundary 
Layer 

Most of the methods pro- 
posed and actually used are 
extensions, frequently very 
direct, of those perfected in 
stationary flow. Two types of /I 5 
methods are distinguished, with- 
out comparison: methods of solu- 

tion of local equations and 
methods of solution of global 
equations . 

3.2.1. Proposed Methods 

In the first type of method, 
the basic equations to be solved 
are local equations of continuity 
and momentum, which are derived 
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from Navler-Stokeo equatlono (1) and (2). The inotantaneous velocity 
is broken down into the average of the set and the turbulent fluc- 
tuation u'^. After applying this breakdown in equations (1) and (2), 
forming the average of the set of equations obtained and using the 
boundary layer hypotheses, the following is obtained 


"bx, 




bx C 






( 36 ) 

( 37 ) 


Different turbulence schemes have been tested. The simplest are 
the "zero equation models," eddy viscosity [^15-48, 7] or the mixing 
length models [10-153 . These models are Identical to those used in 
statlonai'’y flow. For example, the mixing length model is written 


- < 





( 38 ) 


where A/^x,t) is the same function as in stationary flow 

1- 0.025 iS./'iiil l.'i 

S " '0.0‘i^ S ' 

Near the wall, corrections are introduced to allow for viscosity 
effects . 


Several other authors have used a one equation model for the 
Reynolds stress [4, 42, 32] . A two equation model for the kinetic 
energy and dissipation also have been used [10-153 > In these models, 
only the convection term is modified, by including the time derivative. 
For example, the two equation model is written 


3k 

3t 
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4. U,'v'> = Cfj. 
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( 39 ) 

( 40 ) 


In the second type of method, the local equations are integrated 
over y. In general, they use the integral equation of momentum 
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and, frequently, an auxiliary equation, which can be the Integral 
equation of continuity 

Ut 5ot (42) 

or the integral equation of average kinetic energy 




I ii r'’*’ 


e* if 





(43) 


or, again, an equation of momentum, obtained by multiplying the local 
equation of momentum by y 






(44) 


Whatever the auxiliary equation used, it is necessary to have 
available supplementary relationships to solve the integral equations, 
since they bring out an excessive number of unknowns. These closure 
hypotheses are relationships between the integral quantities brought 
out by integration of the local equations. 


McDonald and Shamroth solve equations (4l) and (43). The supple- 
mentary relationships are obtained, by assuming that the average set 
velocity profile can be represented at each instant by a Coles profile. 
The coefficient of dissipation was calculated, by expressing friction 
by means of a mixing length scheme. 


In the method proposed by Kuhn and Nielsen, the basic systems are 
equations (4l) and (44), which are solved after linearization, by 
assuming that the unsteady perturbation is weak. The supplementary re- 
lations were obtained by assuming that, at each instant, the velocity 
profile is described by a Coles profile. 

M.H. Patel has analyzed two methods. Each solves linearized equa- 
tions, obtained by assuming a small unsteady, harmonic perturbation. In 
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the firot method, the average flow Is oalculated by a stationary method, 
for example, the method of Green. With the average flow known, the 
oselllatory perturbation can bo calculated, based on two hypotheses: 

1. at each Instant, the wall friction can be calculated by a stationary 
rule (the Ludwieg-Tillmann law, for example)] 2 , the in phase and out of 
phase components of the velocity are represented by 

Wjj t V 

Patel has obtained these profiles from his experiments. These 
expressions were utilized to calculate the perturbations of the 
Integral densities. R and S thus become the main unknowns in the Integral 
equations, and the solution provides their values. 

In the second method analyzed by Patel, it simply is assumed that, 
at each Instant, the supplementary relations used in the method of 
Green remain valid. 

The method proposed by Cousteix et al [11-133 likewise is an 
extension of the method first established in stationary flow [31] . The 
basic system is made up of the Integral equations of momentum (4l) and 
continuity (42) . The closure relationships are obtained, after deter- 
mination and analysis of the similarity solutions (similar in spirit to 
the stationary solutions of Palkner and Skan of laminar flow) ■ It is 
assumed that the deficit velocities obey a similarity rule of the form, 
(Ug-U)/U^=P' (y/(S (x,t) ) , and the friction is calculated by a mixing 
length model. It then is shown that the resulting family of profiles 
is strictly identical to that determined in stationary flow. It /1 6 

depends uniquely on the Clauser parameter: 

The set of closure relationships necessary to solution of the integral 
equations is then deduced: 

one rule for the ratio 

0 V F, (i|7) 


(45) 

(46) 
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Qne a:’Ulo for Uio wall friction 

4- - 7“ 4| A + .‘DV<5 ) 

one rule for the drive oooffieient 
^ 0 |. (,^ 1 , 

With reopect to the stationary case, only the relationship for 
the drive coefficient is modified. 

D** and P are functions of G, determined by the similarity solu- 
tions 

P= o,o?<(G » 

3.2.2. Application of Different Methods 

Despite the large number of publications on calculation of unsteady 
turbulent boundary layers, it is difficult to compare the performance 
of the various proposed methods. Frequently, the methods have been 
applied to purely theoretical cases, and the various authors have not 
taken up the same cases. L.W. Carr recently proposed standardization of 
these examples. Without being devoid of interest, such exercises, mean- 
while, are of only relative value. One method is compared with another, 
and it is quite clear that no universally valid method exists in tur- 
bulent flow. Only systematic comparison to experiment permits deter- 
mination of the field of validity of a method. Therefore, here, we limit 
outselves to reporting some results of such comparisons. 

At present, the most complete example is that of a boundary layer 
of a flat plate exposed to a harmonic perturbation. A comparison with /I 7 
available experimental data is presented in Pig. 20. We note that the 
Strouhal number is not the unique parameter of the solution. The 
Reynolds number and the amplitude of the fluctuations can play a part. 


(H8) 

(H9) 
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Fis* 20a. Calculation-experiment comparison: average zero 

pressure gradient j phase of 6-j^. 



Fig. 20b. Calculation-experiment comparison: zero 

average pressure gradient i amplitude of 6^. 











lb appoaro bhab, for Inbomodiiabo Sbj:'Ouhal nurabora In bhe 1 b© 5 
i*aneo, bho vanlouD oaleulabiono are nob boo aifforont fnom oxporimonb . 

Poi* Xango SbrouhaX numbeno, bhore io a major* dlffleulby. Ib Id bho 
modeling of burbulonee in bhe vioeoao undor’lying layer*. Aebually, v/hon 
bho Sbrouhal number* Id VGr*y lar*go, bhe unobcady offoebo ar*o ineluded in 
a very bhin layer* near* bho wall, and prediebion of bhom Id eDoonbially 
aDDOoiabed with modeling of the underlying layer, and bhio problem Id 
part of that of the Inberactlon bebv/oen bhe forced puloabion and tur- 
bulence, since the forced frequency Id in bhe range of bho character isblc 
frequencies of the turbulence. Acharya has carried out experiments 
with pulse flow in a duct, in v/hlch the amplitude and phase variations 
of the velocity are located mainly in the underlying layer. The cal- 
culation tests he presents effectively indicate difficulties. 
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Pig. 21. Calculation-experiment comparison: Intensity 

of turbulence. 


We present a comparison of turbulence intensity and turbulent 
pressure profiles we have obtained in Pigs. 21 and 22. With the exper- 
imental dispersion, especially of <u'v'>, taken into account, it is 
difficult to state which of the methods used is better. However, it 
appears that the mixing length model predicts the deformation of the 
<u'v'> profiles more poorly, especially at Instants 9 and 12, located 
at the external velocity maximum. 
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Pig. 22, Galculation-experlmont comparison: turbulent 

proDouro . 


Pigs. 23 and 2k show a comparison of the experiments of M.H. Patel. 

The different calculations were carried out with integral methods. Up /I 8 

to the highest Strouhal number (6.65), good agreement with experiment 
is obtained, the phase deviations certainly are not significant, with 
the difficulty of determining them experimentally with precision taken 
into account. In any case, a small dispersion must be noted in the 
theoretical results for the largest Strouhal numbers. 

Finally, Pig. 25 presents a comparison of the first results we 
obtained with the average pressure gradient . The calculations, carried 
out by the integral method, are interrupted at the station where a /19 

return flow appears. Actually, extension of the calculation area in 
the recirculation region requires that certain conditions at the limits 
along the boundary downstream of the area be taken into account. 

3.2.3. Breakaway problem. Formation of Singularities in the 

Boundary Layer Calculation. 

The breakaway problem raises numerous other questions, and it still 
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Pig. 23. Comparison of experiments of M.H. Patel, 




V; 



Pig. 24. Comparison of experiments of M.H'. PateT. 
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lo not resolved, even in laminar 
flow. 

The only certain proposition 
is negative. Breakaway is not 
necessarily connected to point 
Cf^=0, as in stationary flow. It 
must be emphasized that confusion 
frequently has occurred over the 
word, breakaway, since some authors 
have, by definition, connected this 
word to the fact that Cf=0. 

Rather, it must be reserved to 
signify catastrophe (breakaway, 
breakdown) for the boundary layer 
equations, by which the very no- 
tion of boundary layer ceases to 
be valid. 


Several authors have proposed 
breakaway criteria, principally in 
laminar flow. We first mention 
that of Moore, Rott and Sears [39] j 
established for the case of sta- 
tionary flow with a mobile wall 
Breakaway develops at a point where the velocity Is zero and 
where the gradient also Is zero. This problem Is a special case of un- 
steady flow because. In a reference system connected to the wall, the 
flow Is quite unsteady, and It has been proposed [39] to use a similar 
criterion In the general case. However, numerical studies have shown 
that a singularity does not always occur at this point. 


Pig. 25. Calculation-experiment 
comparison; pulsed turbulent bound- 
ary layer with average pressure 
gradient . 


(Up?^0) . 


A more global approach to this problem, always in laminar flow, 
has been studied by Shen and Nennl, who found the existence of a singu- 
larity, by a noncoincidence condition between the boundary layer and the 
potential flow. This condition is expressed by the fact that the vertical 
velocity becomes not limited at the boundary of the boundary layer. 

^1 



They have shown that the wall friction behavior gives an indication of 
this condition. Besides, they have shown that the wall friction obeys 
a Burgers type equation 


jv - q e 


^ 1 / 
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h^, hg and h^ are functions of x and t. 


Therefore, it clearly appears that point 3?p=0 does not, a priori, 
play any special part, except in stationary flow (3 =0), where a Gold- 
stein type singularity then is found for r^, which^^ varies as the 
square root of x. In unsteady flow, the formation of singularities must 
be found elsewhere, in the form of a discontinuity by development of a 
shock wave. 


Comparable results have been found by Cousteix et al in the tur- 
bulent case, by analyzing the properties of the Integral equations 
(equation of continuity and of momentum), to which closure relation- 
ships 1,^7), (^8) and (49) are added. It leads to the following 
conclusions : 


the system always has two real characteristic directions 
and (X=dx/u^dt)j it is hyperbollcj 

one of the directions, is always such that 0<Xj^<l^ 

the 2nd direction, is positive for values of the form 

factor less than a critical value H (H ^ 2 . 6 ); it is negative for H>H . 

It also is shown that point H=H , in practice, is confused with 
point Cj,=0. Therefore, the result is that the formation of singularities 
is not connected with point 0^=0. When Xg becomes negative, this 
signifies that information is transmitted from downstream to upstream. 

The development of singularities occurs by the formation of shock 
fronts, across which the characteristics of the boundary layer are 
discontinuous. It is clear that the hypotheses themselves of the 
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boundary layer, then, are defective. It can be thought that, aa in 
stationary flow, coupling with the external flow must be resorted to, 
in order to avoid these singularities. 

To support this hypothesis we consider the case of a separated 
boundary layer, calculated as the limit of a transient development, 
when time becomes infinite. To simplify, we analyze the case of a 
theoretical flow in a unidimensional diffuser (Pig. 26). 

For t<0, the external velocity 
if constant, and the boundary layer 
obeys the stationary rules. At t= 

0, the velocity is modified dis- 
continuously . Ug decreases linearly 
in the 0<x<0.6 range, and is 
constant for x>0.6. For t>0, two 
calculations have been carried out. 
In the first, without coupling, the 
forced velocity at t=0 is maintained 
independent of time. Then, when t 
increases, a discontinuity in the 
evolution of 6^ is found, which 
corresponds to the formation of a 
shock front. 

In the second calculation, an 
attempt at coupling between the 
boundary layer and the nonviscous 
fluid was made, simply by using a unidimensional section principle 

^ l_UeCS^2S.)..o (51) 

l)t £>«■ 

where S(x) is the diffuser cross section such that, at time t=0, the 
velocity distribution satisfies Ug(S-26^) =Cte . 

In the second calculation, therefore, the velocity distribution 
develops as a function of time, according to Eq. (51). Under these 
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Pig. 26. Calculation of station- 
ary boundary layer with recircu- 
lation airhole by unsteady method 
with coupling to external flow 
taken into account. 


oon(5itlons, a perfectly stable stationary solution results. It Is of 
Interest to note that this solution includes a return flow offset, and 
that it shows no sign of singularities at points C|,=0 . We stress the 
fact that this could only be obtained, by means of coupling between the 
nonviscous fluid and the boundary layer. 

We point out that Nash and Scruggs have carried out turbulent 
boundary layer calculations, by using a transport equation model. In 
the theoretical examples with v/hlch they deal, they also observe the 
formation of singularities, if the external velocity is forced. By 
using a pressure gradient reduction procedure when the thickening of 
the boundary layer is too great, they show that the singularity can be 
avoided by maintaining a zone with return flow. 

We also point out that Briley and McDonald have proposed a boundary 
layer calculation technique with a breakaway offset, by an unsteady 
method, while taking account of coupling with the nonviscous fluid. 

If. Conclusions 


Recent experiments and those which are underway are providing 
indispensable data on the behavior of oscillating turbulent boundary 
layers, which have been missing so far, and they can be the basis for 
deciding the validity of the calculation methods . 

Meanwhile, some areas still remain partly unexplored. In the quite 
moderate Strouhal range (S<5), the effect on the viscous underlying 
layer of unsteadiness must be better defined. This problem perhaps Is 
not too crucial for Strouhal numbers, since the proposed methods of cal- 
culation give results which, as a whole, are quite consistent with each 
other and with those of experiments, at least when the average pressure 
gradient Is zero. Por higher Strouhal numbers (on the order of 20), the 
problem Is much greater. Actually, unsteady perturbations show up In a 
layer, the thickness of which decreases when the Strouhal number Increases 
and can be limited to the thickness of the underlying layer. Therefore, 

It Is thought that prediction of the coefficient of friction and of Its 
amplitude and phase depends mainly on good modeling of the turbulence In 
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the underlying layer. Under such conditions, a second problem is raised: 
for these Strouhal numbers and the customary Reynolds number, the fre- 
quency Imposed on the flow, within the range of characteristic frequencies 
of the turbulence. Interaction between the harmonic perturbation and 
turbulence can then exist. 

Breakaway and the formation of singularities in the boundary layer 
calculation remain unsolved problems, even in laminar flow. Several 
numerical studies, in laminar, as well as turbulent flow, have shown that 
singularities can develop Inside the boundary layer, but interpretation 
and theoretical study of them are Incomplete. These singularities fre- 
quently appear to clearly oppose the validity of the hypotheses used to 
establish the boundary layer equations. Nevertheless, it can be hoped 
to remove these difficulties in a number of cases, by turning to coupling 
with the external flow. 
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